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Table 1. NGC 1501 Observing Sites 







Pixel 


Read Noise 




Exposure 


Telescope 


Detector 


Size 


(electrons) 


Filter 


00 


Kitt Peak 0.9-m 


Tektronix 512 x 512 


0'.'77 


8 


Stromgren y 


90 


Lick 1-m 


GEC 288 x 385 


0'.'80 


9 


7500 A 


120 


Okayama 0.9-m 


Hamamatsu 1024 x 1024 


0'.'20 


15 


Johnson V 


120 


Wise 1-m 


RCA 320 x 512 


0'.'87 


54 


Johnson V 


120 


Wendelstein 0.8-m 


Thomson 384 x 576 


0'.'47 


20 


Bessell R 


120 



Notes to Table 1. 

Observers were as follows: 
Kitt Peak: Bond 
Lick: Stover 

Okayama: Tamura, Malasan, Yamasaki, Hashimoto, Kambe, Takeuti 
Wise: Leibowitz 

Wendelstein: Soffner, Roth, Mitsch 
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Table 2. Table of observations of NGC 1501 for the 1991 global campaign 



Telescope 


Run Name 




Date 


Start Time 


Run Length 








(UT) 


(UT) 


(hh:mm:ss) 


Okayama 0.9m 


ok91-001 


14 


1NOV 


1 001 

iyy i 


10:11:15 


9:12:56 


Okayama 0.9m 


ok91-002 


15 


Nov 


1 001 

lyyi 


12:30:52 


7:14:52 


Okayama 0.9m 


ok9 1-003 


18 


1NOV 


1 001 

lyy i 


13:56:26 


5:51:11 


Wise 1.0m 


wi91-001 


18 


1NOV 


1 001 

lyy i 


16:14:16 


10:52:39 


Lick 1.0m 


H91-001 


19 


1NOV 


1 001 

iyy i 


6:07:29 


5:29:52 


KPNO 0.9m 


kp91-001 


19 


Nov 


1 001 

lyy i 


7:03:58 


5:54:13 


Wise 1.0m 


wi91-002 


19 


1NOV 


1 001 

lyy i 


16:12:40 


8:24:15 


Lick 1.0m 


H91-002 


20 


Nov 


1 001 

lyyi 


3:04:15 


10:24:00 


KPNO 0.9m 


kp91-002 


20 


1NOV 


1 001 

lyyi 


3:25:54 


9:15:59 


Okayama 0.9m 


ok91-004 


20 


1NOV 


1 001 

lyyi 


10:51:50 


8:42:42 


Wise 1.0m 


wi91-003 


20 


Nov 


1 001 

lyyi 


16:09:46 


10:28:57 


KPNO 0.9m 


kp9 1-003 


21 


1NOV 


1 001 

lyy i 


3:17:32 


9:19:01 


Okayama 0.9m 


ok91-005 


21 


1NOV 


1 001 

iyy i 


9:58:58 


9:27:31 


Wise 1.0m 


wi9 1-004 


21 


Nov 


1 001 

lyyi 


21:12:24 


4:09:43 


KPNO 0.9m 


kp9 1-004 


22 


Nov 


1 001 

lyyi 


2:42:46 


9:51:37 


Lick 1.0m 


U91-003 


22 


1NOV 


1 001 

lyyi 


8:00:35 


5:36:00 


Okayama 0.9m 


ok91-006 


22 


Nov 


1 001 

lyyi 


9:43:32 


0:48:13 


Lick 1.0m 


H91-004 


23 


1NOV 


1 001 

lyy i 


2:40:20 


10:50:00 


KPNO 0.9m 


kp9 1-005 


23 


Nov 


1 001 

tyyi 


3:06:07 


9:33:27 


Okayama 0.9m 


ok91-007 


23 


Nov 


1991 


18:09:26 


1:03:13 


TXT 111 

Wendelst. 0.8m 


we91-001 


23 


Nov 


1991 


19:12:50 


10:27:40 


Lick 1.0m 


H91-005 


24 


Nov 


1991 


2:36:25 


7:24:51 


KPNO 0.9m 


kp91-006 


24 


Nov 


1991 


2:39:44 


10:00:14 


Okayama 0.9m 


ok9 1-008 


24 


Nov 


1991 


9:23:14 


10:38:06 


Wendelst. 0.8m 


we9 1-002 


24 


Nov 


1991 


18:24:59 


8:50:19 


Lick 1.0m 


U91-006 


25 


Nov 


1991 


2:32:52 


10:49:08 


KPNO 0.9m 


kp91-007 


25 


Nov 


1991 


3:08:32 


9:33:40 


Wendelst. 0.8m 


we91-003 


25 


Nov 


1991 


16:16:09 


12:54:58 


Lick 1.0m 


U91-007 


26 


Nov 


1991 


2:28:25 


10:50:00 


KPNO 0.9m 


kp91-008 


26 


Nov 


1991 


2:43:14 


9:57:37 


Wendelst. 0.8m 


we91-004 


26 


Nov 


1991 


17:01:32 


12:29:07 


KPNO 0.9m 


k P 9 1-009 


27 


Nov 


1991 


2:34:00 


10:10:08 



Table 3. NGC1501 and its Comparison Stars 



Star a (J2000) S (J2000) V B-V 



NGC 1501 
Comp 1 
Comp 2 
Comp 3 



4 h 06 m 59!44 
4 06 42.54 
4 07 06.19 
4 06 57.94 



+60°55'14'/6 
+60 57 13.6 
+60 57 18.4 
+60 56 46.2 



14.20 +0.71 

13.05 +0.95 

13.81 +0.89 

15.45 +1.57 
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Table 4. Periodicities of NGC 1501 from 1991 Global Campaign Data. 



Frequency 
(MHz ) ' 


Period 

00 




Amplitude 
(mma) 


rp a 

■*■ max 

00 


comments 


866.281 ±0.048 


1154. 


.36 


±0. 


.06 


5.00 


±0.37 


132 ±42 




855.508 ±0.130 


1168. 


.90 


±0. 


.18 


1.84 


±0.37 


1034 ±115 


possible alias 


798.777 ±0.084 


1251. 


.91 


±0. 


.13 


2.66 


±0.36 


872 ±81 




758.463 ±0.044 


1318. 


.46 


±0. 


.07 


5.05 


±0.36 


693 ±45 




728.366 ±0.119 


1372. 


.94 


±0. 


.22 


1.87 


±0.35 


323 ±126 




698.555 ±0.124 


1431. 


.53 


±0. 


,25 


1.80 


±0.36 


550 ±137 




661.089 ±0.121 


1512. 


.66 


±0. 


.28 


1.77 


±0.36 


967 ±146 




567.947 ±0.071 


1760. 


.73 


±0. 


.22 


3.15 


±0.36 


102 ±96 




528.275 ±0.088 


1892. 


.95 


±0. 


.32 


2.55 


±0.36 


1778 ±128 




500.209 ±0.116 


1999 


.16 


±0. 


.46 


1.92 


±0.36 


198 ±178 




191.029 ±0.059 


5234. 


.81 


±1. 


.62 


3.73 


±0.36 


4896 ±238 





a Timc of first maximum after BJDD 2448571.50367 = Oh UT 1991 Nov 14 
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Table 5. Archival Photometry of NGC 1501 



XYUll 


OLdl L Uaie 


oiaii J- line 


1 In vofi f\ t~i 

UUI dllOIl 




(VI) 


(VI) 


[Il.IIlIIl.boJ 


k P 87-01 


29 Dec 1987 


4:23:46 


2:25:13 


kp87-02 


30 Dec 1987 


4:23:22 


4:49:00 


kp88-01 


2 Jan 1988 


7:00:34 


3:12:10 


k P 88-02 


27 Nov 1988 


3:28:17 


9:36:59 


kp88-03 


28 Nov 1988 


1:26:43 


11:39:40 


kp88-04 


29 Nov 1988 


1:21:24 


11:32:23 


k P 89-01 


8 Oct 1989 


5:50:55 


6:23:22 


kp89-02 


9 Oct 1989 


6:04:07 


6:19:30 


k P 89-03 


15 Oct 1989 


4:59:51 


6:40:32 


kp89-04 


16 Oct 1989 


7:01:36 


5:39:53 


kp89-05 


21 Oct 1989 


8:42:10 


2:24:34 


kp90-01 


1 Jan 1990 


1:46:18 


9:15:00 


kp90-02 


2 Jan 1990 


1:17:52 


10:01:47 


k P 90-03 


29 Nov 1990 


2:18:41 


11:07:24 


k P 90-04 


30 Nov 1990 


1:42:05 


11:36:33 


kp90-05 


3 Dec 1990 


6:24:54 


6:13:56 


kp90-06 


6 Dec 1990 


7:01:34 


5:40:38 


kp90-07 


7 Dec 1990 


6:52:37 


5:51:30 
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Table 6. Dominant Periodicities of NGC 1501 observed in archival runs. 



Runs 


Frequency 


Amplitude 




(/*Hz ) 


(mma) 




:569.14 ±0.22 


7.7 ±1.6 


kp87-01,02,kp88-01 


:667.97 ±0.21 


9.6 ±1.6 




795.30 ±0.19 


9.4 ±1.6 




p" o o /in i a o/t 

533.49 ±0.37 


5.5±0.9 




fpTi AO I /~i I A 

558.03 ±0.40 


5.2±0.9 




656.04 ±0.18 


11.6±0.9 


kp88-02-04 


707.94 ±0.34 


6.2±0.9 




729.46 ±0.30 


6.9±0.9 




864.70 ±0.68 


3.0±0.9 




900.24 ±0.73 


2.8±0.9 




:524.92 ±1.20 


5.0±1.6 




557.96 ±0.55 


10.7±1.6 




657.39 ±0.61 


9.5±1.6 


kp89-01,02 


708.08 ±1.40 


5.3±1.6 




730.40 ±1.30 


5.9±1.6 




:768.5l ±1.90 


3.1±1.6 




833.19 ±1.00 


5.8±1.6 




865.15 ±1.00 


5.8±1.6 




532.46 ±0.72 


8.5±1.8 


k P 89-03,04 


568.12 ±0.68 


9.0±1.8 




666.25 ±0.48 


12.2±1.7 




532.50 ±0.09 


6.2 ±1.1 




558.30 ±0.06 


8.7 ±1.1 


kp89-01-05 


655.78 ±0.05 


9.9 ±1.1 




:708.29 ±0.11 


5.0 ±1.2 




:730.31 ±0.11 


5.0 ±1.2 




:768.93 ±0.12 


4.4 ±1.1 




865.89 ±0.09 


5.1 ±1.1 




:533.84 ±0.78 


5.3 ±1.1 




:557.64 ±0.76 


5.5 ±1.1 


kp90-01,02 


655.78 ±0.45 


9.1 ±1.1 




707.85 ±0.79 


5.2 ±1.2 




728.97 ±0.89 


4.7 ±1.2 




865.87 ±0.62 


6.5 ±1.1 




655.75 ±0.91 


5.1 ±1.3 


k P 90-03,04 


695.86 ±0.47 


9.8 ±1.3 




758.52 ±0.61 


7.4 ±1.5 




864.33 ±0.76 


5.9 ±1.3 
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Table 7. A Model for the Pulsation Spectrum of NGC 1501 a 



Frequency* 3 Period b m An Model Period Pmodei — fobs 



( M Hz) 


00 






00 


00 


500.2 


1999 





22 


2003 


+4 


528.3 
533.3 


1893 




-1 


17 


1891 


-2 


558.0 
(562.9) 
568.0 


(1776) 


+1 

(0) 
-1 


12 


1780 


+4 


656.0 
661.1 
666.3 


1513 








1512 


-1 


698.6 
(703.3) 
708.0 


(1422) 





-4 


1423 


+ 1 


(723.4) 
728.4 


(1382) 


(0) 

-1 


-6 


1379 


-3 


758.5 
(762.5) 


(1312) 


+1 

(0) 


-9 


1312 





798.8 
(803.8) 


(1244) 


+1 

(0) 


-12 


1245 


1 


866.0 


1155 





-16 


1156 


1 



a P = 1512.4s + 22.297s x An. 

b Frequencies and periods for m = modes that are inferred from 
the observed frequency spacings are indicated in parentheses. 
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Asteroseismological Observations of the Central Star of the Planetary 

Nebula NGC 1501 



Of 



Howard E. Bond 1,2 , Steven D. Kawaler 3 , Robin Ciardullo 4 , R. Stover 5 , T. Kuroda 6 , T. Ishida 6 , 
T. Ono 6 , S. Tamura 7,8 , H. Malasan 7,9 , A. Yamasaki 7,10 , O. Hashimoto 7,11 , E. Kambe 7,10 , 
M. Takeuti 7 ' 8 , T. Kato 12 ' 13 , M. Kato 12 ' 13 , J.-S. Chen 14 , E. M. Leibowitz 15 , M. M. Roth 16 ' 17 , 
VO i T. Soffner 16 , 

and 

; W. Mitsch 16 

m : 

: 

ABSTRACT 

CN ■ 

We report on a global CCD time-series photometric campaign to decode the pulsa- 
; tions of the nucleus of the planetary nebula NGC 1501. The WC4 central star is an 

extremely hot, hydrogen-deficient, "O Vl"-type object, with some spectroscopic charac- 
teristics similar to those of the pre-white-dwarf PG 1159—035 stars. NGC 1501 shows 
pulsational brightness variations of a few percent with numerous individual periods rang- 
ing from 19 to 87 minutes. The pulsation amplitudes and periods are highly variable, 
suggesting a complex pulsation spectrum that requires a long unbroken time series to 
resolve. To that end, we obtained CCD photometry of the central star over a two-week 
period in 1991 November, using a network of observatories around the globe. We ob- 
tained nearly continuous coverage over an interval of almost one week in the middle of 
the run. With this data set, we have identified 10 independent pulsation periods, rang- 

>- . ing from 5235 s down to 1154 s. The pulsation modes changed amplitude significantly 

• i— i 

X 
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during the course of the run, indicating either real amplitude variations, or that the 
modes are not fully resolved over the two-week interval. We find strong evidence that 
the modes we see in this star are indeed nonradial g-modes. The ratios of the frequen- 
cies of the largest-amplitude modes agree closely with those expected for modes that are 
trapped by a density discontinuity in the outer layers. This conclusion is strengthened 
by including single-site observations of this star, obtained during previous years, in our 
analysis. We offer a model for the pulsation spectrum that includes a common period 
spacing of 22.30 s and a stellar rotation period of 1.17 days; the period spacing allows 
us to assign a preliminary asteroseismological mass of 0.55 ± O.O3M . However, several 
factors complicate the analysis. Aside from the proximity of the rotational splitting to 1 
cycle per day, this frequency splitting corresponds closely to period spacings near 20 sec- 
onds near the dominant frequencies of the star. Thus, the period spacing and frequency 
spacings are nearly degenerate. 

Subject headings: planetary nebulae: central stars — stars: pulsation — stars: evolution 
- planetary nebulae: individual (NGC 1501) 
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1. Introduction 

A new and powerful tool for probing the depths 
of white dwarfs is asteroseismology. Stars that are 
multimode nonradial pulsators can reveal their in- 
ternal structure through accurate measurements of 
the periods of their photometric (or radial-velocity) 
variations. The prototype of the pulsating hot white 
dwarfs, PG 1159-035 (GW Vir), has yielded a wealth 
of information about its interior, through analysis of 
its light curve as obtained by the Whole Earth Tele- 
scope (WET); see Nather et al. (1990) for a detailed 
description of WET. The detailed analysis of WET 
observations of PG 1159—035 reported by Winget 
et al. (1991, hereafter "WWET") and Kawaler & 
Bradley (1994, hereafter KB) illustrates what such 
stars can teach us; this star shows an unbroken se- 
quence of dozens of pulsation modes that match ex- 
pectations from stellar evolution and pulsation the- 
ory remarkably well. Other GW Vir stars, though 
they show fewer modes than PG 1159—035, have also 
been explored in detail through asteroseismology with 
WET data (i.e., PG 2131+066; Kawaler et al. 1995, 
PG 1707+427; O'Brien et al. 1997). In addition, 
PG 0112+200 shows a rich and decodable pulsation 
spectrum (O'Brien et al. 1996). 

As illuminating as the seismological study of the 
four known pulsating GW Vir-type white dwarfs has 
been, planetary- nebula nuclei (PNNs) that pulsate 
have an equally enormous potential as probes of stel- 
lar evolution. Nine PNNs are now known to be non- 
radial pulsators. The first two to be discovered were 
the nuclei of K 1-16 (Grauer & Bond 1984) and Lo 4 
(Bond & Meakes 1990), and six more were discovered 
in a photometric CCD survey of O VI and PG 1159 
PNNs as summarized by Ciardullo & Bond (1996). 
The total is brought to nine by the discovery of a very 
faint planetary nebula around the pulsating PG 1159- 
type star RX J2117. 1+3412 by Appleton, Kawaler, & 
Eitter (1993). (In addition, the PNN-like O VI star 
Sanduleak 3, although lacking a visible nebula, is a 
pulsating variable; Bond, Ciardullo, & Meakes 1991.) 

Unfortunately, PNNs are more difficult to observe 
than the pulsating PG 1159 stars for several reasons. 
One major difficulty is that the pulsation periods of 
the PNNs range from about 700 s to more than 5000 s, 
and thus slow transparency variations over a given 
night make disentangling the true stellar pulsation 
from variations in the sky transparency difficult. An- 
other problem is that, by definition, PNNs are sur- 



rounded by nebulosity. For photoelectric aperture 
photometry this reduces the measured pulsation am- 
plitude, and causes almost intractable problems due 
to spurious variations introduced by seeing and small 
guiding errors. Still another frustration is that in the 
best-studied pulsating PNN to date, K 1-16, a given 
pulsation frequency apparently changes significantly 
from night to night, and changes amplitude irregu- 
larly. Despite heroic efforts (for example, see Grauer 
et al. 1987), the pulsations in K 1-16 have still not 
been fully resolved. While such instability could re- 
sult from intrinsic variations in the star, it could also 
result from beating between modes that are closely 
spaced in frequency, or from confusion resulting from 
diurnal aliases or overlapping rotational splittings (see 
below). This seems to be a property of the entire 
class; most if not all of the pulsating O VI central 
stars show complex and variable behavior (Ciardullo 
& Bond 1996). 

Most of the observational problems listed above 
are solved by global monitoring campaigns with CCD 
cameras. CCD imaging systems are capable of mea- 
suring the relative magnitudes of PNNs on time 
scales as short as the exposure plus readout time, 
while the two-dimensional format, which can include 
simultaneously-observed comparison stars, avoids the 
problems associated with the nebula and variable sky 
transparency. In addition, by using a network of 
CCD-equipped telescopes around the globe, it is pos- 
sible to obtain nearly continuous data on a specific 
PNN for a week or more, thus greatly reducing the 
ambiguities introduced by the daily gaps in single-site 
data. 

In this paper, we report the results of a global ob- 
serving campaign on the central star of the planetary 
nebula NGC 1501. NGC 150l|] was first discovered 
to be a variable star in 1987 by Bond & Ciardullo 
(1993). It is classified as a WC4 central star, and 
shows prominent emission lines of O VI in the stel- 
lar spectrum (Smith & Aller 1969; Heap 1982; Kaler 
& Shaw 1984). Using the Hell Zanstra temperature, 
Stanghellini et al. (1994) estimate that NGC 1501 has 
an effective temperature of logT c ff = 4.91 + 0.03; from 
its visual magnitude and estimated distance, they ob- 
tain logL/L© = 3.33 ± 0.15. Comparison of its po- 
sition in the H-R diagram with evolutionary models 
(see Figure 9 of Stanghellini et al. 1994) indicates 



In this paper we will refer to central stars of planetary nebulae 
by the names of the nebulae. 
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that the mass of NGC 1501 is w 0.56M©. NGC 1501 
is very hot, rich in C and O, and hydrogen deficient. 
It is similar to the PG 1159 spectral class of pre- 
white-dwarf stars (see, for example, Werner et al. 
1991; Werner 1995), although evidently somewhat 
less evolved. NGC 1501 pulsates on a time scale of 
about 25 minutes, with an amplitude of up to several 
percent. In general, therefore, it appears to possess 
typical properties found in the known pulsating PNNs 
(Bond et al. 1993; Ciardullo & Bond 1996). Because 
of the spectroscopic similarities between these stars 
and the "naked PG 1159 stars" it is presumed that 
there is a direct evolutionary link, in the sense that 
the central stars are the precursors of the PG 1159 
stars. NGC 1501 is thus an early link in the chain of 
evolution between stars with active nuclear burning 
(AGB and post-AGB stars) and stars that are des- 
tined simply to cool and fade (white-dwarf stars). 

NGC 1501 is probably the best northern PNN 
candidate for a global monitoring campaign, since 
it is relatively bright (broadband optical magnitude 
of about 14), and because its northerly declination 
(+61°) makes it observable for nearly 11 hours per 
night from a single northern-hemisphere observatory 
during November. 

We are entering into a very exciting time in astero- 
seismology and in the study of PNNs. The NGC 1501 
results presented in this paper represent the first de- 
tailed analysis of the nonradial pulsations of a PNN. 
In the next section, we describe the organization of 
our CCD network for the NGC 1501 campaign, the 
observing procedures, and the resulting time-series 
data. We then describe the results of our analysis 
of the 1991 campaign data, and the pulsation spec- 
trum of NGC 1501. With the aid of our older archival 
data on NGC 1501, we then construct an initial model 
of the underlying pulsation spectrum, and obtain the 
stellar mass and rotation period. We conclude with 
an attempt to understand this pulsation spectrum in 
terms of models of PNNs and PG 1159 stars. 

2. Observations 

2.1. CCD Photometry 

For 14 days in 1991 November, we monitored the 
central star of NGC 1501 with 0.8- to 1.0-m telescopes 
at five sites around the globe: Kitt Peak National 
Observatory (KPNO) and Lick Observatory in the 
United States, Okayama Astrophysical Observatory 
in Japan, Wise Observatory in Israel, and Mt. Wen- 



dclstein Observatory in Germany. (Four additional 
sites also attempted to observe the object during the 
campaign, but due to weather and technical problems, 
their contributions could not be used.) At each site, 
the object was observed throughout the night via a 
series of short (<~ 2 minute) exposures with a CCD 
detector, and a filter which accepted visual or red 
light. In all cases, the field of view of the CCD was 
large enough to include NGC 1501 itself and three 
nearby comparison stars (shown in Figure 1). A sum- 
mary of the telescopes and observing setups used in 
this campaign is presented in Table 1. Table 2 de- 
tails the active observing time of each observatory. 
Good weather prevailed at many of the sites, so that 
between 1991 November 18 and November 28, useful 
data were obtained for 173 hours out of a possible 216. 
For example, all nine scheduled nights at KPNO pro- 
duced usable data (albeit through varying amounts 
of thin cirrus and haze, which are of little concern for 
differential CCD photometry). It may be noted that 
the campaign was centered around full moon, which 
is not a problem for CCD observations of a star this 
bright, and made it easier to obtain telescope time. 

Our data-reduction techniques are described in de- 
tail by Ciardullo & Bond (1996), but are summarized 
here. Since each night of observing produced be- 
tween 200 and 300 CCD images, a frame-compaction 
algorithm was first used to reduce the data to a 
manageable size and format. After debiasing and 
flat-fielding the data, we extracted 32 x 32-pixcl re- 
gions surrounding NGC 1501 and its three compar- 
ison stars from each individual frame. [Because of 
the small (O'.'2/pixel) plate scale of the Okayama ob- 
servatory data, the frames from this telescope were 
first binned 4x4 before the object extraction took 
place.] The extracted regions were then moved into 
a set of "packed" pictures, each containing data from 
16 original frames. Although this operation eased the 
burden of data handling, it in no way affected the 
subsequent data reductions, since the data from each 
original frame were always analyzed independently. 

Photometric reductions were then accomplished 
using DAOPHOT algorithms (Stetson 1987) inside 
IRAF. Each frame's point-spread function (PSF) 
was first defined using the three prc-sclccted compar- 
ison stars. The relative magnitudes of these stars 
and that of NGC 1501 were then computed using 
the DAOPHOT PSF-fitting algorithm with a fitting 
radius of the order of the image full-width at half- 
maximum (FWHM). (Because NGC 1501's nebula 
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is of relatively low surface brightness, no special im- 
age modeling of the nebula was needed [cf. Ciardullo 
& Bond 1996]). From these data, NGC 1501's rela- 
tive intensity was measured by differencing the raw 
magnitude of NGC 1501 and the flux-weighted mean 
instrumental magnitudes of the comparison stars. 

Because of its favorable declination (+61°), each 
telescope in our network could observe NGC 1501 
for over 10 hours a night. However, during these 
long runs, the airmass of the object changed sub- 
stantially. This change, coupled with the color dif- 
ference between the PNN and our comparison stars 
(see Table 3 for coordinates and our measured mag- 
nitudes and colors), could, in principle, introduce an 
apparent long-term variation in NGC 1501's bright- 
ness, especially in data taken through a broadband 
filter. We therefore checked for this effect by regress- 
ing NGC 1501's computed instrumental magnitude 
against airmass for each observatory's data, and mea- 
suring the slope of the relation. In all cases, the sys- 
tematic error was £ 0.015 mag per airmass, with a 
one-sigma uncertainty of about half this value. Since 
the effect is small, and in any case the time scales of 
primary interest are significantly shorter than the in- 
dividual runs, we did not make any correction for dif- 
ferential extinction between NGC 1501 and the com- 
parison stars. 

A second correction to NGC 1501's instrumental 
magnitude comes from the effects of seeing. By re- 
stricting our photometry to apertures of the order 
of the FWHM and treating each frame identically, 
any error caused by mis-estimating the sky plus neb- 
ula surrounding NGC 1501 should have been min- 
imized. However, under conditions of variable see- 
ing, small systematic errors in the relative magni- 
tudes can still exist. To investigate this possibility, 
we again performed a least squares fit, this time re- 
gressing NGC 1501's differential magnitude against 
the FWHM of the each frame's PSF. Once again, the 
seeing correction proved to be small, typically less 
than 0.01 mag per arcsec. As a result, we did not 
correct for this in the final analysis. 

2.2. Time-Series Analysis 

After obtaining the intensity of NGC 1501 on each 
individual frame, we computed the "modulation in- 
tensity" by subtracting off, and then dividing by, the 
star's overall mean intensity. After this was done, 
individual points that were far away from the rest 
(generally due to cosmic-ray hits on the variable or 



one of the comparison stars) were edited out, and the 
observation times were reduced to Barycentric Julian 
Dynamical Date (BJDD). The data from each site 
were reduced separately, and then combined into a fi- 
nal light curve. We made no effort to weight the data 
for the different apertures of the telescopes, CCD sen- 
sitivities, or other site-dependent effects. 

A subset of the data is shown in Figure 2 in 24-hour 
strips for seven days at the heart of the run. Figure 2 
includes only data from Okayama, KPNO, Wise, and 
Wendelstein observatories. (Because the Lick data 
largely overlap those from KPNO, including these 
data would make the plot more complicated. How- 
ever, the agreement between observatories is quite 
good. In the analysis below we do include the Lick 
data when determining accurate frequencies.) The 
principal variation of about 1200 to 1500 s is readily 
visible in this plot. The plot also shows the tempo- 
ral coverage of our data; the lack of any long vertical 
stripes without data is an indication that there were 
no strictly periodic gaps in this segment of the light 
curve. The window function for this run, obtained 
by sampling a single- frequency sine wave of fixed am- 
plitude at all of the observed sample times, is shown 
in Figure 3. The window function has 1 cycle day" 1 
alias peaks at several percent of the power of the main 
peak. 

We computed the Fourier transform of the reduced 
data set using a brute-force discrete Fourier trans- 
form, and squaring the amplitude of the complex 
transform to produce the power spectrum. This gives 
the total modulation power (mp) as a function of fre- 
quency The square root of the modulation power 
gives the modulation amplitude (ma) . In many cases, 
we find that plots of modulation amplitude (which we 
call "amplitude spectra" ) are more useful than power 
spectra, as the amplitude spectra de-emphasize strong 
peaks. The representation of pulsation power and am- 
plitude in units of mp and ma is discussed in Winget 
et al. (1994). The amplitude spectrum is shown in 
Figure 4 for the entire frequency range of interest; 
for comparison, the window function is plotted again, 
below the strong peak at 866 /xHz. 

While the amplitude spectrum shows relatively 
clean peaks, in some cases the (small) alias side- 
lobes affect the perceived amplitudes and frequen- 
cies of lower-amplitude modes. Essentially, the low- 
amplitude peaks can be shifted away from their true 
frequencies by power that leaks from the larger-amplitude 
peaks. Also, the amplitude spectrum does not give 
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a quantitative estimate of the uncertainty in the fre- 
quency, amplitude, or phase of a given periodicity. We 
therefore used a linear least-squares fitting of the time 
series photometry to determine the precise frequen- 
cies, amplitudes, and phases of many peaks simul- 
taneously. Careful use of least squares allows these 
quantities to be determined accurately for all modes 
included in the fit (i.e., Kawaler et al. 1995; O'Brien 
et al. 1996); the least-squares values are unaffected by 
the sampling pattern if all modes are included in the 
fit. This procedure also provides a formal estimate of 
the uncertainty of these quantities. 

We used the amplitude spectrum of the entire data 
set to provide initial guesses for the 11 clearest peaks 
for use in the initial least-squares fit. The least- 
squares frequencies match the peaks in the ampli- 
tude spectrum within the errors. These frequencies 
are listed in Table 4. The times of maxima are with 
respect to T = BJDD 2448571.50367, which corre- 
sponds to a barycentric date of 1991 November 14.0. 
The only suspicious peak identified in the power spec- 
trum is the peak at 855.5 /iHz, which appears to be 
close to the —1 cycle d _1 alias of the large peak at 
866.3 ^Hz. Because of its relatively large power, and 
the difference in frequency between it and the true 
— 1 cycle d _1 alias, we suspect that it might be a real 
periodicity in the star. Because of our uncertainty, 
however, we exclude it in the following analysis. 

2.3. Archival Data for NGC 1501 

Since only 10 modes were definitively identified in 
the 1991 global campaign, these data alone are insuf- 
ficient for the kind of asteroseismological analysis per- 
formed on the GW Vir stars. However, other pulsat- 
ing white dwarfs show modes that appear and disap- 
pear on relatively long time scales; when new modes 
appear, they appear at frequencies expected given the 
fact that they are (/-modes. That is, the new modes 
usually fit a pattern of equal period spacing and/or 
rotational frequency splitting. In the case of PNNs, 
Grauer et al. (1987) and Ciardullo & Bond (1996) 
show that the modes seen in PNNs can come and go 
on much shorter time scales than in GW Vir stars. 
If these modes are indeed normal modes of oscilla- 
tion, they must occur at the limited set of frequencies 
fixed by the structure of the star. Therefore, analysis 
of data taken at earlier seasons could reveal additional 
frequencies which can be used in the analysis. 

In the hopes of finding such additional modes, we 
reexamined time-series photometry of NGC 1501 ob- 



tained by Ciardullo & Bond (1996) at KPNO between 
1987 and 1990. The individual observing runs are 
listed in Table 5. Data were obtained using the same 
CCD system as in the present study, and were re- 
duced following the same procedures as for the 1991 
global campaign. Because these were all single-site 
runs, the power spectra of these additional data are 
riddled with diurnal aliases. However, in several of 
the seasons, a few large-amplitude modes can be iden- 
tified unambiguously. The frequencies, periods, and 
amplitudes of the dominant modes in those runs are 
listed in Table 6. 

We note that these data sets show abundant ev- 
idence for additional periodicities, but the aliasing 
problem is severe. To disentangle aliases from the true 
frequencies we followed the prewhitening and synthe- 
sis procedures of O'Brien et al. (1996). In some cases 
we were able to make confident frequency determina- 
tions for these runs. In Table 6, the most insecure 
frequency determinations arc indicated with colons 
(:). A few points are readily apparent. First, the 
power present in the dominant modes from these runs 
is much larger than in the 1991 campaign. Moreover, 
the period with largest amplitude changes from year 
to year; only in the 1991 data is the 866 fiHz mode 
among the largest. Finally, it is clear that NGC 1501 
displayed a number of modes in the recent past that 
were not present during the 1991 campaign. These 
additional modes can be used to constrain models of 
the underlying pulsation spectrum of the star. 

3. Frequency and Period Patterns in Nonra- 
dial g-modes 

NGC 1501's pulsation frequencies, if caused by nor- 
mal nonradial g-modes, should show some systemat- 
ics. It is these systematics which allow us to place 
seismological constraints on the global parameters of 
the star. This section briefly describes some of the 
systematics of nonradial g-modes that are present in 
pulsating white dwarfs and that we might expect in 
the power spectrum of NGC 1501. 

3.1. Frequency Splitting 

One of the observed properties of nonradial oscil- 
lations is that equal frequency spacings are often seen 
when "bands" of power are resolved. For a perfectly 
spherically symmetric star, the frequencies of all non- 
radial modes with a given degree £ and order n are 
identical. However, departures from spherical sym- 
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mctry (which can be caused by rotation, magnetic 
fields, or some combination thereof) can lift the de- 
generacy. The most common example is rotation, 
which splits the degenerate modes into 21 + 1 com- 
ponents (if all values of the azimuthal index m are 
present) as follows: 

V n lm = V n lQ-m- 5 —{l-C r d), (1) 
-frot 

where for white dwarf stars, to within a few percent, 

Cnl = W+~i) (2) 

(Brickhill 1975; WWET 1991). Thus, if NGC 1501 
rotates, we can expect to see triplets split by one half 
of the rotation frequency for I = 1 modes. We follow 
the convention that the lower-frequency component 
of a triplet is identified as the m = +1 component, 
and the higher-frequency component is labeled as the 
m = — 1 component. Such equally spaced (or nearly 
equally spaced) triplets are commonly seen in pul- 
sating PG 1159 stars and in the cooler DB and DA 
white-dwarf pulsators. 

3.2. Period Spacings 

In the limit of high overtones (that is, when the 
number of radial nodes is much greater than £) g- 
modes are approximately equally spaced in period for 
a given value of I and consecutive values of n, i.e., 

n "' = 7fe ( " + e) - (3) 

where e is a small number. The period spacing (IIo) is 
determined largely by the mass of the star, and is only 
weakly dependent on the stellar luminosity and sur- 
face compositional stratification. KB, and references 
therein, discuss how the value of this period spacing 
is used to constrain the mass of GW Vir stars, and 
how departures from uniform spacing allow the deter- 
mination of subsurface structure. Since the pattern 
of equally spaced (in period) multiplets has been seen 
in several of the GW Vir stars with period spacings 
of approximately 20-23 seconds, we might expect a 
similar spacing in the periods of NGC 1501. 

The relatively small number of observed modes in 
NGC 1501, coupled with the long periods and result- 
ing lack of period resolution at the frequency resolu- 
tion set by the run length, makes finding period spac- 
ings a difficult task using the 1991 data alone. As we 
show in the next section, the archival data allow us 
to identify a possible period spacing. 



3.3. Strange Frequency Ratios 

WWET and KB found that the period spacings 
seen in PG 1159 itself were not precisely equally 
spaced. They noted that this could result from mode 
trapping; that is, a subsurface composition transi- 
tion region can cause a density discontinuity which 
acts as a reflecting boundary between the surface and 
the deep stellar interior. KB identified this layer as 
the transition zone between the helium-rich surface 
layer and the carbon/oxygen core. In their anal- 
ysis, KB computed so-called trapping coefficients, 
which are related to the frequency ratio of trapped 
modes. If trapped modes are identified in a star, 
then the period ratio between the modes can be com- 
pared with these coefficient ratios. In fact, the domi- 
nant peaks in the amplitude spectrum of the pulsat- 
ing central star RX J2117. 1+3412 show period ratios 
nearly exactly equal to y/3/2 (Vauclair et al. 1996), 
which is compatible with the coefficients computed 
by KB. The precision with which the period ratios in 
RX J2117.1+3412 match V3/2 is remarkable, much 
better than 1%. This suggests that dominant peaks 
in RX J2117. 1+3412 are indeed trapped modes. 

It is immediately apparent that the same ratio ex- 
ists in NGC 1501. The frequency ratio between the 
peaks at 866 ^Hz and 758 ^Hz is 0.8755, or 1.1% 
larger than y/3/2. The next two peaks in the sequence 
would be at about 656 ^Hz and 569 ^Hz. We do see a 
small (but still significant) peak at 661 /xHz (at a ratio 
0.7% smaller than y/3/2) and a large peak at 568 /iHz 
(at a ratio 0.8% bigger than \/3/2). These observed 
ratios are all much closer to to \/3/2 than seen in 
the theoretical models of KB. Apparently the pulsat- 
ing central stars NGC 1501 and RX J2117.1+3412 
are better calculators of \/3/2 than are the computer 
models! 

It is very interesting to see that a large-amplitude 
mode exists in past data at 656-657 /xHz. Figure 5 
shows the spectrum of the 1988 data, along with a 
representative window function. This frequency falls 
exactly in the \/3/2 sequence ratio; the next fre- 
quency expected in the sequence below 758 fiHz is 

656.8 /uHz, and the next frequency expected above 

567.9 ^Hz is 655.8 fiHz. Thus the archival data 
strongly support the identification of a series of pe- 
riodicities with frequencies near the "magic" ratio of 
V3/2. 

4. Decoding the Pulsation Spectrum of NGC 1501 
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4.1. Identification of Rotational Splitting 

The 656 /iHz periodicity that dominated much of 
the archival data is not present in the 1991 data, but 
there is a modest peak at 661.1 /iHz. The frequency 
difference here is about 5.4 /kHz. In the 1987 data 
and the second half of the 1989 data sets, there are 
clear peaks at 668.0 /xHz and 666.3 /iHz, respectively. 
With the 661 /iHz peak from the 1991 campaign, these 
three peaks, seen in different years, may all be compo- 
nents of a uniformly split triplet with 5v = 10.2 /iHz. 
Triplet structures like this are seen in many of the pul- 
sating GW Vir stars; as described in the previous sec- 
tion, it is the expected signature of rotational splitting 
of t — 1 modes. Armed with evidence of a frequency 
triplet centered near 661 /zHz, we used the archival 
data and campaign data to reveal other evidence for 
frequency splittings of approximately 10 /iHz. 

While the 1991 data show an unaliased peak at 
567.9 /tHz, the strongest peak in the runs kp89-01 
and kp89-02 is clearly at 558.0 /iHz. This same fre- 
quency is present in 1988 data. Figure 6 shows this 
region of the power spectrum of the 1988 runs. To 
try to decompose this portion of the spectrum in 
the presence of 1 cycle/day aliases, we follow the 
technique of power-spectrum synthesis described in 
O'Brien ct al. (1996). This procedure allows us to 
find the minimum number of independent frequencies 
that best reproduce the observed power spectrum. In 
the top two panels of Figure 6, we show the observed 
power spectrum as a thick line. In these panels, a thin 
line shows a synthetic power spectrum created using 
two sinusoids with amplitudes, phases, and frequen- 
cies obtained by least-squares fits to the light curves. 
The two frequencies used in each synthesis are marked 
in the figure. The rest of the peaks in the panel arise 
solely from the window patterns associated with the 
input frequencies and gaps in the data; no noise was 
added in this simulation. In the topmost panel of 
Figure 6, the higher- frequency peak has a frequency 
close to that seen in the 1991 data (w 568 /tHz). 
The fitting and synthesis procedure reveals that the 
two frequencies do not adequately reproduce the ob- 
served power spectrum, indicating that either one or 
both frequencies are incorrect, or that more frequen- 
cies are present. However, if we fit the light curve by 
setting the higher frequency at ss 558 /iHz (second 
panel), the entire region is very well matched. Based 
on this agreement, we are confident that our iden- 
tification of a 558 /iHz periodicity in the 1988 data 
is correct. Thus there is a pair of peaks separated 



by about 10 /iHz centered at 562 /iHz. This is (to 
within the frequency errors) entirely consistent with 
the frequency spacings in the 656-666 /iHz region. 

The archival data support two additional frequency 
bands with the identified splitting near 530 /iHz and 
700 /iHz. The 1991 data show a peak at 528.3 /iHz; 
the archival data show a peak at 533.5 /iHz (in 1988; 
cf. Figure 6), and at 532.5 /iHz (in 1989). Thus, 
we have evidence for a pulsation mode approximately 
5 /iHz higher than in 1991. Though there is no evi- 
dence of a third peak in this frequency range, the data 
are consistent with a frequency splitting of modes 
with the same £ and n but consecutive values of m. In 
1988, 1989, and 1990 a strong peak was also evident at 
a frequency of 708 /tHz; this peak lies almost exactly 
10 /iHz above the 698.6 /iHz peak seen in 1991. 

The mean value of the rotational splitting between 
the m = +1 and m = — 1 modes in NGC 1501 
is 9.9 /iHz ; equations (1) and (2) then imply that 
NGC 1501 rotates with a period of 1.17 days. If this 
9.9 /iHz splitting is truly present in NGC 1501, then 
it would explain why this is such a difficult star to 
analyze from a single site! The 1 cycle/day alias sep- 
aration is 11.6 /iHz; data from a single site would 
need to cover at least 11 days to begin to resolve the 
true power from the alias. Since the modes present in 
NGC 1501 appear to come and go on this same time 
scale, identification of the true peaks from a single 
site is nearly impossible. 

4.2. Period Spacings 

While the number of modes is still a bit small, we 
can use the measured frequencies in the archival runs, 
along with the 1991 global campaign, to try to model 
all of the observed frequencies in terms of modes that 
are equally spaced in period and rotationally split in 
frequency with a total splitting of w 10 /iHz. Un- 
fortunately, to search for period spacings, we need to 
identify the m = components of the triplets in those 
cases where all three modes are not present. Compli- 
cating matters further, the period spacings are not 
expected to be precisely uniform, as mode trapping 
by a subsurface composition transition region could 
cause small variations in the period spacing with pe- 
riod, as in the pulsating PG 1159 stars (KB). 

To compute the period spacings, we begin by deter- 
mining the periods of the m = modes. The triplet 
near 656 /iHz tells us that the m = central frequency 
can be estimated when we see two modes separated 
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by 10 /iHz. From the frequencies in Tables 4 and 6, we 
see that m = periods are available for the 563 /iHz 
group, the 661 /iHz group, and the 703 /iHz group. 
The differences in these periods are 263 s and 91 s; 
these should share a common multiple which itself is 
a multiple of the intrinsic period spacing. The ratio 
of these is 2.89; that is close enough that we will call 
it 3; thus the period spacing is either approximately 
90 seconds or some integer factor (i.e., 45, 30, 22.5, 
18, or 15 seconds). 

To reduce this degeneracy, we can appeal to the 
theoretical calculations of periods and period spacings 
in PG 1159 stars by KB. Assuming that these arc I = 
1 modes, and that NGC 1501 is a normal planetary- 
nebula central star, we can eliminate large values for 
the period spacing. Values above 25 s for the period 
spacing would imply a mass M < 0.5M Q , well below 
the minimum mass needed for a single C/O white- 
dwarf star. 

We next consider the mode at 758.5 /iHz seen in 
the 1991 data. With no other nearby frequency seen 
in any of the archival data we cannot claim an iden- 
tification of m. If it is m = 0, then the period is 104 
seconds shorter than for the m = component of the 
703 /iHz group. This is not commensurate with the 
other two period intervals. If this mode is actually 
m = — 1, then the m = component would have a 
period 94.9 s longer than the 703 /iHz group — which 
is also not commensurate. However, if the m = 
component is 5 /iHz higher in frequency (i.e., at a fre- 
quency of 763.5 /iHz), then the m = period for this 
mode is 110 s shorter than for the 703 /iHz group. 
This number is compatible with a period spacing of 
22 s. 

Since an integer number of spacings must separate 
the periods of the modes, the difference between 1776 
seconds and 1312 seconds corresponds to 21 inter- 
vals of 22.1 s. The subintervals between successive 
m = modes are 12 intervals of 21.9 s, four intervals 
of 22.8 s, and five intervals of 22.0 s. With a re- 
fined estimate of the period spacing interval, we can 
continue this analysis and attempt to identify values 
of m and the number of intervals between successive 
modes (An) for each of the frequencies seen in the 
1991 data set (and several from the archival runs as 
well). Finally, a least-squares fit of the periods and 
An's allows identification of the period spacing. This 
spacing is 22.30 seconds. Table 7 presents this model 
for the pulsation spectrum of NGC 1501. In all cases, 
the model matches the observed periods to within 4 s. 



Variations this large (and larger) are present in the 
pulsating PG 1159 stars that have been studied by 
WET, so the model appears to be satisfactory. 

Using equation (5) of KB, and noting that this 
equation contains a sign error on the exponent of the 
luminosity term, we can use the 22.30 s period spacing 
to estimate the mass of NGC 1501 as 0.55 ± 0.03M Q . 
The quoted uncertainty reflects the fact that the lu- 
minosity of NGC 1501 is uncertain; while Stanghellini 
et al. (1994) estimate an uncertainty of only about 
35%, we allow that the luminosity could be anywhere 
between w 5000L Q and w 100i Q to retain a conser- 
vative error estimate in our mass determination. In 
addition, KB derived the relationship between mass 
and period spacing for PG 1159 stars, which are more 
highly evolved than NGC 1501; the period spacing 
does tend towards an asymptotic value at high lumi- 
nosities in their Figure 2, however. 

5. Summary and Conclusions 

With a coordinated global observing campaign, we 
have obtained a high-resolution power spectrum of 
the pulsations of the central star of the planetary neb- 
ula NGC 1501. The star showed 10 clear periodicities 
during this observing run, with periods ranging from 
5200 s down to 1154 s. The largest-amplitude pul- 
sations occur between 1154 s and 2000 s. Additional 
data from prior years reveal that NGC 1501 changes 
its pulsation properties dramatically over month to 
year time scales. Past observations have found it pul- 
sating with amplitudes twice as large as seen in 1991. 
The pulsations present in prior years in some cases 
occur at the same frequencies as in 1991, but in other 
cases the pulsations appear at different (but not ran- 
dom) frequencies. We note that as a PNN, NGC 1501 
must be evolving relatively rapidly; with many modes 
seen at the same frequencies in subsequent years, we 
are optimistic that further analysis of the data may 
yield a measurement of the evolutionary change in 
one or more of the stable modes. 

Most of the pulsation frequencies are accounted 
for by a simple model which includes just two input 
parameters. We attribute the 9.9 /iHz spacing ex- 
hibited by several multiplets to rotational splitting, 
implying that the star's rotation period is 1.17 days. 
Moreover, there is a nearly uniform underlying period 
spacing of 22.30 s between the various multiplets, im- 
plying an asteroseismological mass for NGC 1501 of 
0.55 ± 0.03M Q . This mass is in close agreement with 
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the mass determined from its position in the H-R dia- 
gram and comparison with evolutionary models (i.e., 
Stanghcllini ct al. 1994). Despite its unusual spec- 
trum, then, NGC 1501 is a garden-variety PNN in 
terms of its mass; the mass we find is similar to other 
estimates of PNN masses. It also is consistent with 
the star being a progenitor of the PG 1159 spectro- 
scopic and pulsation class. 

Our success in analyzing NGC 1501 is muted by 
the significant observational difficulties presented by 
this star. NGC 1501 changes its pulsation behavior on 
time scales that prevent detailed monitoring. A global 
campaign is required to ensure that the 10 fiRz rota- 
tional splitting is resolved from diurnal aliases, but 
maintaining such a vigil over the months necessary 
to look for changes in pulsation frequencies and am- 
plitudes would be an enormous task. Moreover, in 
addition to the alias problem, there is an unfortunate 
degeneracy in the pulsation spectrum of this star. A 
period spacing of 22.30 seconds corresponds to a fre- 
quency difference of 10 fiHz, when the periods are 
near 1500 s (i.e., at frequencies near 666 fiHz). Thus 
the rotationally split low-frequency components of the 
power spectrum overlap with the consecutive overtone 
spacing for periods of 1500 s and longer. This may 
in part account for the difficulty faced in resolving 
this star, but unlike diurnal aliasing (which also in- 
troduces similar frequency splitting), observing from 
around the world does not help. 

NGC 1501 is but one of nine PNNs that are now 
known to show nonradial (/-mode pulsations (Ciar- 
dullo & Bond 1996). Data have already been obtained 
for RX J2117.1+3412 by the WET collaboration, and 
are currently being analyzed. We also have data on 
Sanduleak 3 that are currently being examined. With 
this paper's demonstration that NGC 1501 is under- 
standable in terms of well-known principles of astero- 
seismology, it is our hope that the similarities between 
it and the other members of the class will accelerate 
analysis and interpretation of their pulsations. The 
seismological study of this fleeting stage of stellar evo- 
lution has now begun. 
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Fig. 1. — Sample full-frame image of NGC 1501 
and three comparison stars, obtained through a 
Stromgren y filter with the KPNO 0.9-m telescope. 
North is at the top, east is on the left, and the image 
is approximately 4.'3 wide. 

Fig. 2. — Combined light curve of NGC 1501 for seven 
days during the middle of the 1991 global campaign. 
Each horizontal panel represents 24 hours of data, 
with the UT date indicated to the right of the panel. 
Data from KPNO, Okayama, Wise, and Wendelstein 
are included in this plot; since the data from Lick 
overlap the KPNO data almost exactly, those data 
points were not plotted to retain clarity. Data from 
about 3 to 13 hours UT are from KPNO; from 10 to 20 
hours are from Okayama; from 17 to 2 hours are from 
Wise; and from 20 to 5 hours are from Wendelstein. 

Fig. 3. — The spectral window of the global cam- 
paign, computed by sampling a noise-free single- 
frequency sinusoid at the same times as present in 
the data. The top panel shows the power spectrum of 
the window function, while the bottom panel shows 
the amplitude spectrum. 

Fig. 4. — The power spectrum and amplitude spec- 
trum of the combined data from the global campaign. 
The top panels show the power spectrum (and the 
window function) while the bottom panels show the 
amplitude spectrum (and the corresponding window 
function) . 

Fig. 5. — The power and amplitude spectra for the 
central star of NGC 1501, as observed at Kitt Peak 
on three consecutive nights in 1988 November. The 
spectra and windows are as in Figure 4. Note that 
the dominant peak lies at 656 /iKz. 

Fig. 6. — Enlargement of the power spectrum of the 
1988 data from Kitt Peak. The solid line in each panel 
is the data, and the thin lines correspond to synthetic 
power spectra with the indicated input frequencies. 
See text for details of the synthesis technique. The 
second panel shows a good match to the observations. 
The bottom panel represents the window function for 
these runs. 
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